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The origin of strong magnetic fields in the Universe can be explained by amplifying weak seed fields 
via turbulent motions on small spatial scales and subsequently transporting the magnetic energy 
to larger scales. This process is known as the turbulent dynamo and depends on the properties of 
turbulence, i.e. on the hydrodynamical Reynolds number and the compressibility of the gas, and 
on the magnetic diffusivity. While we know the growth rate of the magnetic energy in the linear 
regime, the saturation level, i.e. the ratio of magnetic energy to turbulent kinetic energy that can be 
reached, is not known from analytical calculations. In this paper we present the first scale-dependent 
saturation model based on an effective turbulent resistivity which is determined by the turnover 
timescale of turbulent eddies and the magnetic energy density. The magnetic resistivity increases 
compared to the Spitzer value and the effective scale on which the magnetic energy spectrum is at its 
maximum moves to larger spatial scales. This process ends when the peak reaches a characteristic 
wavenumber fe* which is determined by the critical magnetic Reynolds number. The saturation 
level of the dynamo also depends on the type of turbulence and differs for the limits of large and 
small magnetic Prandtl numbers Pm. With our model we find saturation levels between 43.8% and 
1.3% for Pm ;§> 1 and between 2.43% and 0.135% for Pm ^ 1, where the higher values refer to 
incompressible turbulence and the lower ones to highly compressible turbulence. 


I. INTRODUCTION 

Magnetic fields play an important role in shaping the 
Universe. For understanding the formation of stars and 
galaxies it is thus crucial to know the strength of mag¬ 
netic fields, their distribution in space, and their evolu¬ 
tion in time. Observations in the local Universe indicate 
that magnetic helds are strong. In fact, the energy den¬ 
sity of magnetic fields is often comparable to the thermal 
one and the one of cosmic rays, at least on spatial scales 
above 1 kpc [I] - a phenomenon known as energy equipar- 
tition. Moreover, magnetic fields are observed over a 
huge range of scales from planets [2] and stars [5], to 
interstellar clouds [4], galaxies [5], and potentially the 
intergalactic medium [5]. 
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Generation of magnetic fields is possible already in the 
early Universe [7]. During inflation tiny magnetic fluc¬ 
tuations could expand into large-scale magnetic fields. 
However, with flux freezing in the cosmic expansion the 
field strength decreases rapidly leading to very weak seed 
fields [5] . Likewise other generation mechanisms produce 
extremely weak magnetic fields, as for example held gen¬ 
eration in cosmological phase transitions [SlIID]. Opposed 
to the cosmological origin of seed helds, they can also be 
generated in plasma processes like batteries. For instance 
the Biermann battery m can produce helds of the or¬ 
der of 10“^® — 10“^^ G on a comoving scale of 10 kpc 
m which is roughly 20 orders of magnitude below the 
10“® G helds observed in present-day galaxies m- A 
recent study by Schlickeiser [14] suggests that aperiodic 
plasma huctuations can result in magnetic helds of the 
order of 10“^® G in the intergalactic medium and 10“^° 
G in a protogalaxy. Thus, magnetic seed helds need to 
be amplihed very efficiently during the evolution of the 
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Universe. 

Let us assume that the intergalactic medium has a field 
strength of Bi = 10“^° G and a number density of 
ni = 10“® cm“^. Under the condition of flux freez¬ 
ing, the field strength increases only to a value of B 2 = 
( 711 / 712 )“^/^ Bi = 10“^® G in spherical gravitational col¬ 
lapse during the formation of a galaxy with a number 
density of 712 = 1 cm“^. Pure compression of the field 
lines during the collapse of a halo is thus insufficient to 
amplify these weak seed fields. Another possible ampli¬ 
fication mechanism is a large-scale galactic dynamo that 
converts the kinetic energy from the galactic rotation into 
magnetic energy and simultaneously orders the fields m- 
However, this process operates on long timescales compa¬ 
rable to the one of galactic rotation. In order to explain 
the observed magnetic fields in galaxies, the large-scale 
dynamo would require stronger seed fields than predicted 
by the generation mechanisms discussed above. 

There is probably only one universal mechanism by which 
the tiny seed fields can be amplified to the observed 
present-day galactic values: the small-scale or turbulent 
dynamo. This type of dynamo converts kinetic energy 
from turbulent motions into magnetic energy by ran¬ 
domly stretching, twisting, and folding the field lines. 
As the turbulent dynamo operates initially on very small 
spatial scales, in fact on the diffusive scales, it is asso¬ 
ciated with very short timescales. In the kinematic dy¬ 
namo phase the magnetic energy grows exponentially at 
large growth rates Tk, increasing the magnetic energy on 
small spatial scales by many orders of magnitude. Subra- 
manian m is one of the first authors who discusses the 
turbulent dynamo as a crucial ingredient for generating 
strong magnetic fields in galaxies. He argues that the 
latter needs to operate to provide a sufficiently strong 
seed field for a large-scale galactic dynamo which in turn 
provides an explanation for the origin of the observed 
fields correlated over kpc length scales. Recent semiana- 
lytical calculations as well as numerical simulations 
[IS] suggest indeed that the turbulent dynamo operates 
efficiently already in young galaxies leading to unordered 
fields with a strength of the order of 10“® G. 

It has been shown that the turbulent dynamo operates 
under a large range of different physical conditions as 
long as the magnetic Reynolds number 


exceeds a critical value Rmcrit- Here V is the velocity on 
the forcing scale L of turbulence and rj is the magnetic 
diffusivity. Dynamo amplification is possible from very 
small [IMSSl to very large [24H26| ^ magnetic Prandtl 


^ We refer here to the preprint paper Subramanian m as it in¬ 
cludes a derivation of the Kazantsev theory for the turbulent 
dynamo as well as technical details on the calculation of the 
kinematic growth rate. 


numbers 


1C Rm 

( 2 ) 

Pm = - = -—, 
r] Re 

where ic is the viscosity and 


VL 

( 3 ) 

Re=- 

V 

is the hydrodynamic Reynolds number. 

Pm describes the 

separation between the viscous and the resistive scale. 

C = Re-i/P+’’) L 

( 4 ) 

and 


= Rm-PP+’’) L, 

( 5 ) 


and can be written as Pm = with d be¬ 

ing the slope of the turbulence spectrum in the inertial 
range. The amplification process is fundamentally dif¬ 
ferent for the two extreme cases of the Prandtl number. 
For Pm ^ 1 the dynamo operates fastest on ty leading 
to a growth rate of Tk oc [26]. Amplifica¬ 

tion for the case of Pm ^ 1 is most efficient on with a 
growth rate of Tk oc [53]. Note that there 

is a strong dependence on the turbulence spectrum via 
its slope d which differs for different types of turbulence. 
While small-scale dynamo action is most efficient in in¬ 
compressible Kolmogorov turbulence with 7? = 1/3 [37] . 
it has been demonstrated that the dynamo can also oper¬ 
ate in supersonic turbulence [ 531|331 [ 531123121 ] j i-e. in the 
regime of high Mach numbers, where the gas is compress¬ 
ible and "d can reach values up to 1/2 for Burgers-type 
turbulence |5D] . 

As soon as the fields have significant back reactions on 
the fluid via the Lorentz force, the peak of the magnetic 
energy spectrum is shifted to larger spatial scales in the 
non-linear dynamo phase. This modification of the spec¬ 
trum can be explained by an increase of the magnetic dif¬ 
fusivity. Subramanian m models an effective diffusivity 
which depends on the magnetic energy density and a re¬ 
sponse timescale. With the assumption that the latter is 
the eddy turnover time we can calculate the diffusion rate 
on a given spatial scale. Dynamo amplification comes to 
an end when the diffusion rate equals the growth rate. 
The magnetic energy continues to increase on larger spa¬ 
tial scales via stretching, twisting, and folding by larger 
turbulent eddies shifting the peak of the spectrum [32) . 
Saturation of the dynamo occurs when the spectral peak 
reaches a maximum spatial scale which is in our model 
determined by an intrinsic property of dynamo action, 
the critical magnetic Reynolds number Rmcrit- We note 
that an alternative approach has been suggested in litera¬ 
ture [33] that models saturation as a result of the velocity 
statistics becoming anisotropic with respect to the local 
magnetic field. 

The fraction of kinetic energy that is converted into 
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magnetic energy by the turbulent dynamo, the satura¬ 
tion level, has been estimated from magnetohydro dy¬ 
namical (MHD) numerical simulations. For incompress¬ 
ible Kolmogorov-like turbulence, the saturation levels are 
typically 40 % 0[3S]. In the large Mach number regime, 
where the gas is highly compressible, Federrath et al. [29] 
find considerably lower values between 6 x 10“^ % and 
1 % depending on the numerical parameters of the simu¬ 
lations. In this paper we aim to compare these and other 
numerical results from the literature with a phenomeno¬ 
logical model for the saturation of the turbulent dynamo. 
The paper is organized as follows: After a brief review 
on the dynamo in its different phases, the kinematic and 
the non-linear phase, in section |ll| we introduce our phe¬ 
nomenological model of the saturation process in section 
|III| We model the behavior of the magnetic diffusivity 
r] which changes to an effective diffusivity r]eS in the 
presence of a strong magnetic field and we discuss the 
consequences of this for the resistive scale and the 
effective magnetic Reynolds number Rmes = VL/rjeS. 
Dynamo amplification shows fundamental differences for 
the regime of large and small magnetic Prandtl num¬ 
bers. Thus, we discuss the two extreme cases separately 
in sections III B and III C| In section jlllD we analyse the 
saturation of the turbulent dynamo for typical parame¬ 
ters of numerical simulations in order to make a direct 
comparison. We draw our conclusion in section EYI 


II. AMPLIFICATION OF RANDOM 
MAGNETIC FIELDS 


The turbulent dynamo amplifies magnetic fields by 
randomly stretching, twisting, and folding the field lines 
in turbulent motions under the condition of flux freez¬ 
ing. This mechanism is most efficient on the smallest 
spatial scales of the system which have the shortest eddy 
turnover timescales. 

For large magnetic Prandtl numbers Pm, amplification is 
fastest on the viscous scale given in Q, where turbu¬ 
lent eddies are the smallest and the fastest. In fact, the 
magnetic energy can even be transported below the vis¬ 
cous scale of the system and accumulates on the resistive 
scale £ri given in ([^. This mechanism is illustrated in 
figure [l| While the turbulent eddies stretch, twist, and 
fold the magnetic field on a scale the field lines are 
brought closer together on spatial scales £ < £i, within 
multiple stretching processes. This way magnetic energy 
is transported down to the resistive scale £ri below which 
it dissipates. 

For the contrary extreme of small magnetic Prandtl num¬ 
bers, Pm ^ 1, the resistive scale is above the viscous 
scale. In this case the dynamo operates fastest on £rj 
and is restricted by the turnover rate of the turbulent 
eddies on that scale. With £^ the typical dynamo 
timescales are larger for small Prandtl numbers and we 
expect the dynamo to be less efficient. 




FIG. 1. An illustration of the amplification of magnetic fields 
and the resulting transport of magnetic energy to subviscous 
scales. Motions of the fluid on the viscous scale stretch 
and bend the magnetic field lines to scales I < Iv This way 
anti-parallel field lines are brought closer together and the 
magnetic energy eventually accumulates on the resistive scale 
£^. 


A. Kinematic Dynamo Phase 

The time evolution of the magnetic field is governed 
by the induction equation, 

f)Fi 

= V X (-Utot X Btot) + (6) 

which relates the total field strength Btot to the velocity 
field t^tot and the magnetic diffusivity rj. By separat¬ 
ing the mean from the turbulent components, i.e. t>tot = 
(v) + V and Btot = {B) + B, one can derive evolution 
equations for the mean-field and the turbulent dynamo. 
An analytical treatment of the turbulent dynamo is 
possible within the framework of the Kazantsev theory 
[T5l [36] . A central assumption of this model is that tur¬ 
bulence is incompressible, where the velocity fluctuations 
V scale with the eddy size £ as w oc However, in 

astrophysical environments turbulence is usually driven 
by high-Mach number flows leading to a different scaling 
[SMI- In this work we use 

V ex. £^ (X k~^ ( 7 ) 

with £ = 27r/fc, where k is the wavenumbeij^ and the free 
parameter £} ranging from 1/3 [57] up to 1/2 for highly 
compressible turbulence [ 301115 ]. Being interested in the 
statistical properties of the turbulent dynamo, we model 
the two-point correlation function {vi{ri,t)vj{r 2 , s)) of 
the fluctuating velocity field for a general type of turbu¬ 
lence which enters via the variable £). The spatial posi¬ 
tions ri and r 2 are related by £ = |r 2 — ri|. If the corre¬ 
lation function of the magnetic field {Bi{ri,t)Bj{r 2 ,t)) 
is rewritten with a separation ansatz proportional to 
^(r) exp(2Ft), the growth rate F is governed by the 


^ We note that we will analyse the spectra in the following mainly 
in Fourier space. Here we use the word “scale” as well for 
wavenumbers k as for spatial lengthscales i. 
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Kazantsev equation [5B] . 

- = -FkV'W- (8) 

d r 

The functions Kdiff(?’) and U{r) depend on 
{vi{ri,t)vj{r 2 ,s)) and 77 . With a model for the 
velocity correlation function one finds generalized ex¬ 
pressions for Kdiff(?') and U{r). The Kazantsev equation 
can be solved analytically in the limits of large and small 
magnetic Prandtl numbers (§ with the result of [221 US] 

f 7k.Pm»i y , Pm » 1, 

V (9) 

[ 7 k,Pm<i ^ R-m^^ ’?)/(!+’?) ,Pni<l, 

with the proportionality constants 


eddy timescale of the current peak scale [35] . The growth 
rate in this so-called non-linear phase is given as 

where we determine the free parameter 7 ni self consis¬ 
tently in the following. The evolution of the magnetic 
field strength on the forcing scale is now independent of 
the Reynolds numbers, but still depends on the slope of 
the turbulence spectrum d. Schleicher et al. [13] find that 
the field strength on the turbulent forcing scale evolves 
proportional to 


III. THE FLOW OF ENERGY OVER SPATIAL 
SCALES AND SATURATION OF THE 
TURBULENT DYNAMO 


3047? -h 163 

7k,Pm»l = gQ (10) 

for large magnetic Prandtl numbers, and 

^ exp - 1) - 2 ) 

\ 91'^ J 

with gi = 56 — 103'!? and g 2 = (79 — 157'!?)'!? — 25 for 
small Prandtl numbers. The numerical values of these 
proportionality constants are 7k,Pm»i = 4-41 — 5.25 and 
7k,Pm<i = 0.0268 — 0.00543 for the slopes of the turbu¬ 
lence spectrum 7 ? = 1/3—1/2. These solutions agree well 
with numerical solutions of equation Q [33] . 

We note that various assumptions are made in the 
Kazantsev theory [36]. One simplification is that the 
turbulent velocity field is delta-correlated in time. Bhat 
and Subramanian m drop this assumption and derive 
a generalized Kazantsev equation for a finite correlation 
time Tcorr- The resulting growth rates are reduced by 
a factor 45/56 Tcorr- In the classical Kazantsev model 
there is neither kinetic nor magnetic helicity included 
which would appear as an additional term in the cor¬ 
relation functions. Malyshkin and Boldyrev [48] study 
the influence of kinetic helicity on small-scale dynamo 
amplification and find that it does not affect the fastest 
growing bound eigenmodes by much for incompressible 
turbulence. 


B. Non-linear Dynamo Evolution 

With the exponential growth, the fluctuating magnetic 
field quickly becomes strong enough to influence the ve¬ 
locity field via the Lorentz force. Then amplification on 
the dissipative scale comes to an end and the magnetic 
energy is shifted to larger spatial scales, typically on the 


At some point the strength of the magnetic back reac¬ 
tions on the velocity field via the Lorentz force is compa¬ 
rable to the dynamo amplification. When this is the case 
on all scales of the system, the magnetic field amplifica¬ 
tion comes to an end, i.e. the dynamo is saturated. 

After we present a model for the magnetic diffusivity 
when the magnetic field approaches saturation on a given 
wavenumber fc, we discuss the evolution of the magnetic 
energy during the different amplification stages and what 
determines their end. In each phase we compare the am¬ 
plification rate with the dissipation rate, and finally de¬ 
rive the magnetic energy at saturation. 


A. Modification of the Magnetic Diffusivity for 
Strong Fields 


Subramanian m suggests a model for calculating the 
saturation energy density of the magnetic field based on 
describing the change of the velocity field. In fact, he 
introduces an effective magnetic diffusivity, 

Stt 

7eff = 7 + 2a-^77i, (13) 


where 77 is the usual microscopic resistivity, i.e. for exam¬ 
ple the Spitzer resistivity, and the parameter a = t/{A^ p) 
with T being the response time of the system and p be¬ 
ing the fluid density. The second term in the effective 
diffusivity (13) is proportional to the magnetic energy 


density m. When the dynamo reaches saturation the ad¬ 
ditional drift component dominates over the microscopic 
resistivity 77 and we obtain from equation (13) 


^eff 


4 r 

--m. 

3p 


(14) 


The timescale t on which the fluid reacts to the mag¬ 
netic field generated by the dynamo should be compara¬ 
ble to the timescale of the turbulent eddies. As the eddy 
turnover time is different on different lengthscales, also 
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the saturation process should be scale dependent. We 
model here the response time on a scale k as 


T{k) = 


27T 

kvk 


27T 

vi4 


k 




(15) 


where we used equation 0 to find the dependence on 
the forcing scale and velocity, /cl and V. 

The rate at which the magnetic energy is dissipated can 
be estimated via the dissipative term in the induction 
equation, 


— = Ves{k)y^B, (16) 


where we approximate d/dt « Tdis and « (27rfc ^) 
The dissipation rate on a scale k is then given by 


^dis{k) 


r]es{k) 

(27rfc“^)2 


(17) 


B. The Limiting Case of Large Magnetic Prandtl 
Numbers 


1. Kinematic Phase for Pm S> 1 


In the kinematic dynamo phase the turbulent energy 
spectrum E{k) is, per definition, not affected by the mag¬ 
netic field. The total initial turbulent kinetic energy is 
thus determined by 

Co = / E{k) dk 
JkL 

= (l - , (18) 

where we used v{k]^) = V and v(k^) = V [k^/k^Y (see 
equation]^ and equation Q to write this in terms of the 
Reynolds number. The total magnetic energy density in 
the kinematic phase rn-k can then be determined from 
the integral over the magnetic energy spectrum in the 
kinematic phase Mk(fc), 


mk 



Mk(A:) d/c. 


(19) 


A Fourier analysis of the Kazantsev theory [SOI leads to 
the following form of the magnetic spectrum in the kine¬ 
matic phase, 

MYk) = a^,k^/‘^Ko (k/kr,), ( 20 ) 


where ak is a constant and KYk/kY) is the modified 
Bessel function of the second kind. The slope which 
characterizes the spectrum at small wavenumbers, 
is known as the Kazantsev slope. Numerical simulations 
have confirmed this shape of the magnetic energy 
spectrum Bhat and Subramanian m show 


that the 3/2 slope of the magnetic spectrum remains 
even if the turbulence has a finite correlation time. 
In the kinematic phase the magnetic energy increases 
exponentially in time with a peak on the resistive scale, 
until: 


Tk — rfiis(77lk, ^r/)- 


( 21 ) 


A comparison of the rates yields the total magnetic en¬ 
ergy at the end of the kinematic phase 


mk = - 7k,Pm»l 




pV^. 


( 22 ) 


For deriving equation (22) we have employed the growth 
rate ([^ and the dissipation rate 0 , where we inserted 
(141 at uik and . Further we used the definition of 
the viscous scale 0 and the relation Pm = 

With (221 we determine the proportionality constant ak 


in the spectrum (20) in the kinematic phase. The re¬ 
sulting spectrum Mk{k) is illustrated in figure]^ by the 
curves (a) and (b) for a hydrodynamic Reynolds num¬ 
ber of Re = 10® and a magnetic Prandtl number of 
Pm = 10^. We present the two extreme cases of tur¬ 
bulence with the Kolmogorov type (i? = 1/3) in the up¬ 
per panel and the Burgers type {d = 1/2) in the lower 
panel. The curve (a) represents a spectrum within the 
kinematic phase, while curve (b) shows the spectrum at 
the end of the kinematic phase. Integration of curve (b) 
over k yields the magnetic energy (22). 


2. Non-Linear Phase for Pm ^ 1 


In the non-linear dynamo phase the magnetic energy 
continues to increase at rate (12), while the peak of the 


spectrum moves to larger spatial scales, i.e. to smaller 
wavenumbers k. It is important to note that the amplifi¬ 
cation rate is now a function of the scale k on which the 
dynamo operates the fastest. This scale is initially the 
viscous scale fcj^, i.e. the scale on which the turbulent eddy 
timescale is the shortest. With the growth of magnetic 
energy, the magnetic diffusivity (13) increases, shifting 


the effective resistive scale on which the magnetic 
spectrum peaks from to smaller wavenumbers. Once 
the peak scale of the spectrum moves to k < k^, amplifi¬ 
cation is fastest on the current resistive scale below 
which the magnetic energy dissipates. 

We thus need to model the non-linear phase in two steps. 
In the first step the effective resistive wavenumber is 
larger than the viscous wavenumber {krj^eff > k^) and the 
growth rate is constant and largest on k^. Once the peak 
wavenumber becomes smaller than the viscous wavenum¬ 
ber {krj^eS < kY amplification is fastest on 

a. First Non-Linear Phase In order to have a steady 
increase of magnetic energy, we demand that the energy 
at the beginning of the first non-linear phase mn\,i is 
equal to the one at the end of the kinematic phase mk. 
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With the field amplification taking place at ki, and dissi¬ 
pation taking place essentially on we find 


r„i(fci,) = rdis(mk, krf) (23) 

at the beginning of the first non-linear phase. The condi¬ 
tion ( |23[) fixes the free parameter in the non-linear growth 
rate (|12|) to 


7nl = 7k,Pm»l- (24) 

Now the peak wavelength krj^eS moves to smaller fc, while 
the current magnetic energy mnii(fcr)eff) is determined 

by 


hnl)^!/) — rdis(mnln(^77,eff); ^r/,eff)- (^b) 


Solving (25) for the magnetic energy and inserting 
constant (24) results in 


4 

(26) 

Note that the magnetic energy increases with the peak 
scale moving to smaller wavenumbers k. The first 

non-linear phase comes to an end when reaches the 
viscous scale k^- At this time the magnetic energy is 


m„u(A:.) = ^ 7k.Pm»i pV^. (27) 



0 1 2 3 4 5 6 7 


We assume that the shape of the magnetic energy spec¬ 
trum is similar to the one in the kinematic phase (201. 
With the peak of the spectrum being at the current re¬ 
sistive scale kri^eS and the normalization changing as a 
function of we model the magnetic spectrum as 

follows: 


Afnl.l (^; ^77,eff a) Q^nl,! (^?7,eff) ^ ^ A^O (^/^?7,eff) • 

(28) 


The normalization constant Q;ni.i(fc?).eff) can be found us¬ 
ing ([^ and 


nTnl.l (^7;,eff) — / Afnl.l (^5 ^?^,efT) dA:. (29) 

JkL 


The resulting energy spectra are illustrated by the light- 
blue lines in figure]^ The curve (d) shows the spectrum 
in the first non-linear dynamo phase at a later point in 
time than curve (c). One can clearly see how the peak 
of the spectrum moves to smaller k in time and how the 
total magnetic energy increases from (c) to (d). 

b. Second Non-Linear Phase Once the effective re¬ 
sistive wavenumber kri^eS moves below the viscous scale 
ku, magnetic field amplification becomes fastest on the 
current The magnetic energy TOni ,2 in this phase 

is then determined by the condition 


logio (fc [fci]) 

FIG. 2. The spectra of turbulent kinetic and magnetic en¬ 
ergy, E{k) and M{k). The input parameters for this plot are 
the hydrodynamic Reynolds number Re = 10®, the magnetic 
Prandtl number Pm = 10® and the slopes of the turbulence 
spectrum = 1/3 and = 1/2 describing Kolmogorov and 
Burgers turbulence. The green dashed line indicates the (ini¬ 
tial) kinetic energy spectrum E{k). The gray lines, (a) and 
(b), show the magnetic energy spectrum M\^{k) at different 
times in the kinematic phase. The two light blue lines, (c) and 
(d), give the spectrum in the first nonlinear phase Mni,i{k), 
while the two dark blue lines, (e) and (f), give the one in the 
second nonlinear phase Mni, 2 {k). The purple line, (g), rep¬ 
resents the magnetic spectrum at saturation Mni, 2 (fc*). In¬ 
dicated as vertical dotted lines are the forcing scale kL, the 
peak scale at saturation A:*, the viscous scale A:„ and the initial 
resistive scale kn- 


which leads to 


tnnl,2(A:r^,eff) 


3 

^ 7k,Pm»l 


kr 


2i? 


pV^. 


(31) 


The peak of the spectrum, 


rnl(^77,eff) — rdis('^nl,2(^77,eff): 


f^nl,2 ) — Q^nl,2 ^ ^ -^0 i (32) 
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continues to move to smaller wavenumbers with the nor¬ 
malization constant Qfni. 2 (fc 7 ;,eff) being determined by 


^nl,2) — / -^17nl,2dfc, (^1^) 

JkL 


when inserting the expression (311 on the left-hand side. 
The dark-blue lines, curves (e) and (f), in figure]^ show 
two representative spectra in the second non-linear dy¬ 
namo phase. While curve (e) shows a spectrum at the 
very beginning of the second non-linear phase, i.e. at a 
point in time when the peak has just moved below ky, 
the spectrum (f) is already closer to saturation. 


3. Saturation for Pm 1 


Saturation occurs when the peak of the magnetic en¬ 
ergy reaches a certain scale /c* which is determined by the 
critical magnetic Reynolds number Rmcrit- R has been 
shown that a minimum Rm is needed for a turbulent 
dynamo to operate which can be interpreted as a mini¬ 
mum separation between the turbulent forcing scale and 
the dissipative scales. This critical magnetic Reynolds 
number can be found by solving the Kazantsev equation 
([^ with r = 0. A typical value for the critical mag¬ 
netic Reynolds number is Rmcrit ~ 100 for Kolmogorov 
turbulence and increases with increasing compressibility. 
For Burgers turbulence Rmcrit ~ 2700 which we use for 
our modeling in this work. A fit formula for the critical 
magnetic Reynolds number is [26j 

Rmcrit(i?) = 88 (tan(2.7 -f 0.2) - 1). (34) 


We note, however, that the values of Rmcrit we use here 
are assumptions based on our previous work. There are 
additional numerical calculations |46| that yield a larger 
value of Rmcrit ~ 32000 for highly compressible turbu¬ 
lence. 

With the additional drift as proposed by Subramanian 
m, an effective magnetic Reynolds number can be de¬ 
fined as 


Rnioff 


VL 

Tj + 16/37r a m 


(35) 


for a given magnetic energy m. Close to saturation Rmotf 
becomes 


Ro^cff (^r/,eff ) 


_ VL _ 

16/3 TT a mnl, 2 (fc» 7 ,eff) ’ 


(36) 


where we can insert the magnetic energy in the second 
non-linear phase (31). It is obvious from equation (36) 
that Rmcff decreases with increasing magnetic energy 
Wni ,2 or equivalently with decreasing The decrease 

of Rmcff comes to an end as soon as it reaches Rmcrit, 
which is the case when reaches a certain scale fc*: 


Rmcfr(A:*) = Rmcrit- (37) 


With equation (36) and the magnetic energy in the sec¬ 
ond non-linear phase (31), we find the following value for 

h- 

k* = (7k.Pm»l Rmcrit)^^^’’"'’^^ fci- (38) 

The numerical values of fc* are listed in table |l] for dif¬ 
ferent types of turbulence. We chose the following rep¬ 
resentative turbulence slopes from literature: The two 
extreme cases, incompressible Kolomogorov turbulence 
with f) = 1/3 [27] and highly compressible Burgers turbu¬ 
lence with -d = 1/2 |30|. If intermittency of Kolmogorov 
turbulence is taken into account the slope of the spec¬ 
trum is i? = 0.35 [52]. In observations of molecular clouds 
typically higher values of the order of d = 0.38 m and 
L = 0.47 [ 35 ] are found. Further we list the slopes re¬ 
sulting from numerical simulations in driven supersonic 
MHD turbulence (i? = 0.37) [ 53 ] and from solenoidally 
{■d = 0.43) and compressively (r? = 0.47) driven super¬ 
sonic turbulence m. 

Using equation (31) the magnetic energy at saturation 
can be calculated as 


TOsat = mnl,2{k*) 

= ^ 7k.Pm>l (7k.Pm»lRmcrit)"^’’^^’’'^^VR^- 

(39) 

The spectral distribution of magnetic energy at satura¬ 
tion is given as 

A7sat(fc) = A7nl,2(fc, ^*) 

= ani, 2 {k*) k^/^ Ko{k/k.), (40) 

where again Kq denotes the Bessel function of second 
kind. We plot the magnetic energy spectra at saturation 
for the two extreme cases of turbulence, d = 1/3 and 
d = 1/2, in figure!^ One major difference between the 
saturation spectra (g) for the different types of turbu¬ 
lence is the value of fc*, which is with a value of 588 fc/, 
much larger for f} = 1/2 than it is for d = 1/3, where 
fc* = lOlkL- Also clearly visible from figure]^ is that the 
total magnetic energy at saturation is much larger for 
Kolmogorov turbulence than for Burgers turbulence. 

We define the saturation level of the turbulent dynamo 
as the ratio of magnetic energy at saturation rusat to the 
kinetic energy at saturation. As the turbulent kinetic 
energy is concentrated on large spatial scales, i.e. within 
the motions of the largest turbulent eddies, the kinetic 
energy spectrum and thus the total kinetic energy is not 
affected much by the growth of the magnetic field. We 
expect the kinetic energy at saturation to be between the 
two extremes of eo and eo — Wsat, leading to the following 
limiting expressions for the saturation level: 

^sat 3 r T3 \—2t?/(i?-|-l) 

- = X 7k.Pm»l (7k.Pm>lRmcrit) (41) 

eo z 

and 

msat _ 3/2 7k,Pm»l (7k,Pm>lRmcrit)~^'^'^^’^''~^^ 

eo ~ jeisat 1 — 3/2 7k,Pm»l (7k,Pm>lRmcrit) ^ 

(42) 














i9 

Rmcrit 

A:* 

Pm > 1 

Pm ^ 1 

?^sat/t0 

Pm > 1 Pm < 1 

^sat/(c0 

Pm > 1 

- msat) 

Pm < 1 

1/3 

« 107 

101 kL 

2.21 kL 

0.304 

0.0238 

0.438 

0.0243 

0.35 

« 118 

104 kL 

2.34 kL 

0.261 

0.0221 

0.352 

0.0226 

0.37 

« 137 

110 kL 

2.51 kL 

0.212 

0.0196 

0.269 

0.0199 

0.38 

« 149 

114 kL 

2.60 kL 

0.190 

0.0182 

0.235 

0.0185 

0.43 

« 227 

135 kL 

2.78 kL 

0.108 

0.0118 

0.121 

0.0120 

0.47 

« 697 

260 kL 

4.15 kL 

0.0410 

0.00457 

0.0433 

0.00460 

1/2 

« 2718 

588 kL 

6.02 kL 

0.0134 

0.00135 

0.0136 

0.00136 


TABLE I. Listed are the slopes of the turbulence spectrum for different types of turbulence (see the discussion in the text) 
and the corresponding critical magnetic Reynolds numbers Rmcrit from [26]. For the different turbulence models we present 
our results for the peak scale of the spectrum fc* at saturation, and the energy ratios msat/eo and msat/(eo — uisat). Here rrisat 
is the magnetic energy density at saturation and eo the initial turbulent kinetic energy. For the values of the kinetic energy, 
we assume that we are in the limit of large hydrodynamical Reynolds numbers, where eo = l/2pV^ is constant. 





model are listed in table |T] for different types of turbu¬ 
lence. Moreover, we illustrate the dependence of the sat¬ 
uration level on the slope of the turbulence spectrum 
1 ? in figure Therefore we insert the fit for the criti¬ 
cal magnetic Reynolds number (34) and shade the area 
between the functions (41) and (42). The upper pur¬ 
ple curve yields the resulting dynamo saturation level for 
large magnetic Prandtl numbers which decreases from 
approximately 30% — 44% at i? = 1/3 to roughly 1.3% at 
d = 1/2. Due to the small values of rusat compared to 
Co at very compressive turbulence, i.e. for r? —>■ 1/2, the 
ratios msa,t/£o and rnsat/(eo ~ Wsat) become comparable 
and the band in figure]^ gets very narrow. 


C. The Limiting Case of Small Magnetic Prandtl 
Numbers 


FIG. 3. The saturation level of the turbulent dynamo as a 
function of the slope of the turbulence spectrum i9. The upper 
boundaries of the bands are given by msat/(eo — nisat), while 
the lower boundaries are determined by msat/co- The upper 
purple band represents the limit of large magnetic Prandtl 
num bers Pm with the saturation levels given in equations 
(411 and (42l. The small Pm limit is illustrated by the lower 
gray band, where the boundaries are equations (511 and (52 1 . 


The saturation level given in equation (41) is the ratio 
of the magnetic energy at saturation rusat to the initial 
turbulent kinetic energy cq. As seen in MHD simulations 
[53|, the kinetic energy spectrum does not change signif¬ 
icantly in presence of a turbulent dynamo, leaving the 
total kinetic energy approximately constant. This can 
be achieved when the dissipation process in MHD turbu¬ 
lence differs from the pure hydrodynamical case and the 
kinetic energy is only partially converted into heat, while 
the remainder is transformed into magnetic energy. If the 
viscous dissipation is similar to the one in hydrodynam¬ 
ical turbulence, energy conservation requires the kinetic 
energy at dynamo saturation to be cq — msat, which leads 
to the saturation level given in equation (421. 

The two energy ratios at saturation gained from our 


For small magnetic Prandtl numbers the resistive 
wavenumber krj is smaller than the viscous one k^,. Mag¬ 
netic field amplification thus initially proceeds fastest on 
kjj as for wavenumbers k > krj magnetic energy is dis¬ 
sipated. When the field becomes strong enough the ef¬ 
fective resistive scale moves from to smaller k 

until it reaches fc* and the dynamo is saturated. We dis¬ 
cuss the evolution of the magnetic energy and its spec¬ 
tral distribution in the following. While the calculation 
is in principle similar to the case of Pm ^ 1 we point 
out the main differences which occur mostly in the non¬ 
linear phase. For Pm ^ 1 the latter is not split into two 
different stages as k^, > kri^eS during the entire dynamo 
amplification. 


1. Kinematic Phase for Pm 1 

The magnetic energy on the resistive scale krj increases 
exponentially at rate ^ until 

Pk = rdis(mk,fc,,). (43) 

Note that the growth rate scales with the magnetic 


























9 
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logic (k [ki]) 

FIG. 4. The spectra of turbulent kinetic and magnetic energy, 
E{k) and M{k). The input parameters for this plot are a 
hydrodynamic Reynolds number of Re = 10® and a magnetic 
Prandtl number of Pm = 10“^. This figure is conceptually 
similar to figure 


Reynolds number Rm for Pm <C 1, while for the opposite 
case of Pm 1 it scales with Re. The condition ( [4^ is 
fulfilled when the magnetic energy reaches 

mk = ^7k.Pm«i pV\ (44) 

The spectral energy distribution can be determined in 
the same way as described in section [III B[ 


.S. Non-Linear Phase for Pm <C 1 

The non-linear phase for the case Pm ^ 1 differs from 
the case of Pm 1. For small magnetic Prandtl num¬ 
bers the amplification proceeds initially on the resistive 
scale. In the transition to the non-linear phase the effec¬ 
tive resistive scale shifts to smaller wavenumbers fc, with 
amplification being always fastest on There is no 

separation into two different non-linear phases. 


With the condition that the magnetic energy at the end 
of the kinematic phase is equal to the initial energy in 
the non-linear phase, we determine the proportionality 
constant 7 ni in the growth rate ( |l^ . With 

rnl(^r7) — ruis(mk, (4^) 

we find 


7nl = 7k.Pm<l- 


(46) 


We note that in equation (451 we are comparing the 


growth rate Tni on with the dissipation rate Pdis on 
kri, while for large Pm we were comparing the growth 
rate on k^, with the dissipation rate on (see equation 

For an arbitrary resistive scale the magnetic en¬ 

ergy in the non-linear dynamo phase is determined by 
the condition 


rnl(^r;,eff) — Tdis(^nl(^77,eff)j ^77,eff)■ 


(47) 


Solution of equation (47) yields the magnetic energy 

2i9 


^nl(^?7,eff) — * 'Tk,Pm'Cl 


^77,eff 


pV^. 


(48) 


This result for the magnetic energy is similar to the one 
for Pm ^ I given in (31) except for the different propor¬ 


tionality constant 7k,Pm<i- 


3. Saturation for Pm 1 


Saturation sets in when the effective resistive scale 
kjj^eS reaches the saturation scale fc*. The latter is de¬ 
termined via the critical magnetic Reynolds number as 
described in section IIIB and is given in equation (381, 


where one needs to replace the constant 7k,Pm>i by 
7k.Pm<i: 

k* = (7k,Pm<i Rmcrit)^^^’’"^^^ kL- (49) 

The magnetic energy at saturation is then determined by 

3 /. \2i9 

mni(A:*) = - 7k,pm<i pV^. (50) 

Again the only difference between the saturation energy 
here and the case of Pm (391 is the proportionality 
constant in the growth rate ([^ which is different for small 
and large Prandtl numbers, see equations (II) and (10). 


For completeness we also give the saturation levels for 
the case of Pm <C 1: 

n^sat 3 r T3 2i?/(i?-|-l) /r-i\ 

— o 7k,Pm<g;l (qk.Pm^lRcUcrit) (bl) 

Co 2 

which is the ratio of magnetic energy at saturation rusat 
to initial kinetic energy eg and 

TOsat _ 3/2 7k,Pm<l (bk.Pm^lRmcrit)”^’’'^^’’"'’^^ 


™sat 1 — 3/2 7k^Pm<l (7k,Pm<lR-IIlcrit) 

( 52 ) 
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FIG. 5. The saturation level of the turbulent dynamo as a 
function of the Mach number M. The crosses are the results 
of the numerical simulations from Federrath et al. for 
solenoidally driven (purple crosses) and compressively driven 
turbulence (gray crosses). The hydrodynamical Reynolds 
number of these simulations is Re ~ 1500 and the magnetic 
Prandtl number is Pm « 2. The black dashed-dotted line at 
low Mach numbers indicates the saturation level in simula¬ 
tions of Haugen et al. [33]. For comparison we include the 
predictions from our model by the blue shaded regions within 
the two limiting cases msat/(£o — f^^sat) and msat/eo (see equa¬ 
tions 41 and 42 1 for different slopes of the turbulence spectrum 
= 1/3 and = 1/2. 


We expect the dynamo to saturate at a value between 
and ( [ 5 ^ . The evolution of the magnetic energy 

f ectra for the case of Pm ^ 1 is presented in hgure 
The exemplary values chosen here are Re = 10® and 
Pm = 10“® and we show the two extreme cases d = 1/3 
(upper panel) and r? = 1/2 (lower panel). The curves (a) 
and (b) show the spectra in the kinematic phase, while 
(b) represents the end of the kinematic phase. Curves 
(c) and (d) are spectra in the non-linear phase and (e) 
is the spectrum at saturation. In comparison to the case 
of large magnetic Prandtl numbers, one clearly sees that 
the magnetic energy spectra are, also at saturation, much 
smaller. In fact, the magnetic spectra are always below 
the initial kinetic energy spectrum for Pm ^ 1. 


D. Comparison with Numerical Simulations 

The turbulent dynamo has been studied intensively in nu¬ 
merical simulations [inilllllMlESlESlES]- In these simu¬ 
lations the full set of MHD equations is solved while tur¬ 
bulence is initiated by a forcing term in the Navier-Stokes 
equation. Different types of turbulence can be simulated 
by applying different types of forcing, i.e. solenoidal forc¬ 
ing, compressive forcing or a mixture of both, and by 
using different Mach numbers A4 which are defined as 
the ratio of the velocity dispersion to the sound speed. 


This way a broad range of different slopes of the turbu¬ 
lence spectrum r? can be studied. 

The dynamo saturation levels found in the high- 
resolution simulations of Federrath et al. |35] are pre¬ 
sented in figure The characteristic numbers of these 
simulations are Re « 1500 and Pm « 2. The ratio of 
magnetic energy to kinetic energy at saturation is plotted 
as a function of the Mach number Ai for solenoidal and 
compressive forcing. The sub-sonic regime of solenoidally 
driven turbulence is comparable to the divergence free 
Kolmogorov turbulence {'d = 1/3), while the supersonic 
regime is dominated by shocks and is better described by 
Burgers turbulence (d = 1/20 We include in figure 
our predictions for the saturation levels using equations 
(411 and ( [4^ as upper and lower boundaries. Note, that 
for = 1 /2the two extremes cannot be distinguished by 
eye in figure The agreement between the predictions 
of our analytical model and the numerical simulations is 
very good for the case of Kolmogorov turbulence, where 
the saturation level is of the order of 40%. In the super¬ 
sonic regime our prediction lies with a value of roughly 
0.64% in between the saturation levels of solenoidally and 
compressively driven turbulence. Further numerical sim¬ 
ulations based on a different code have been performed 
by Haugen et al. [34] who study non-helical turbulence 
without imposed large-scale fields. They use low Mach 
numbers to analyse the dynamo in incompressible turbu¬ 
lence and find a saturation level of roughly 40 % which 
is in agreement with the numerical results for the satu¬ 
ration levels from Federrath et al. [35] and our analytical 
model for the case of incompressible turbulence. We in¬ 
dicate the saturation level found by Haugen et al. [33] by 
the black dashed-dotted line in figure]^ 

A common limitation of simulations is a relatively low 
separation of the dissipative scales, i.e. ki, and k^j, which 
is reflected in the magnetic Prandtl number being of the 
order of Pm « 1. In recent high-resolution simulations, 
Federrath et al. [21] were able to analyse MHD turbu¬ 
lence for a larger range of magnetic Prandtl numbers 
from Pm « 0.1 up to Pm « 10. With Mach numbers up 
to 11, these calculations cover supersonic, highly com¬ 
pressible turbulence. We present the numerically found 
saturation levels as a function of Pm in figure The 
simulations predict an increase of the saturation mag¬ 
netic energy with the Prandtl number. 

We aim to compare these results with the analytical 
model for the saturation of the turbulent dynamo de¬ 
veloped in the previous section. Therefore we determine 
the saturation levels msat/(co) and msat/(eo ~ Wgat) for 
the two extreme cases of Pm 1 and Pm ^ 1. We 


® We note that for the case of f? = 1/2 and a Reynolds number of 
Re = 1500 the saturation scale fc* lies actually above k,^. This 
means that, despites being in the large Prandtl number regime, 
we do only have one non-linear phase. Here saturation takes 
place when the peak of the magnetic spectrum is still at k > . 

This complication does not occur for t) = 1/3 and a Reynolds 
number of Re = 1500. 
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FIG. 6. The saturation level of the turbulent dynamo as 
a function of the magnetic Prandtl number Pm. The pur¬ 
ple crosses are the results of the numerical simulations from 
Federrath et al. |21]. In the analytical model the satura¬ 
tion level is independent of Pm in the limits Pm <C 1 and 
Pm ;§> 1. We shade the region between the two limiting cases 
m.aat/(eo — uisat) and rrisat/eo for different slopes of the tur¬ 
bulence spectrum = 0.4 (dark blue dotted lines), i? = 0.45 
(light blue dashed lines) and = 0.5 (gray solid lines). 


apply a hydrodynamic Reynolds number of Re = 1600 
as given in Federrath et al. [29] and test three different 
slopes of the turbulence spectrum: = 0.4, = 0.45, 

and = 0.5. We expect this range of slopes for super¬ 
sonic turbulence from numerical simulations with high 
Mach numbers |29l HU [54] . The resulting saturation lev¬ 
els are indicated in figure For large magnetic Prandtl 
numbers we find the best agreement between the numer¬ 
ical results and our analytical model for a spectral slope 
of d = 0.45 typical for compressive turbulence. 

In the contrary limit. Pm ^ 1 and even for the regime 
Pm « 1, Federrath et al. |5H| find a decrease of the satu¬ 
ration level with decreasing Pm. In our analytical study, 
on the other hand, we obtain the limiting case of very 
small and very large Pm, without explicitly determining 
the behavior inbetween the two regimes. It is, however, 
clear that the saturation has to drop towards lower val¬ 
ues of Pm in order to fulfill these limits. The discrepancy 
from the simulations at Pm < 1 may arise from the fact 
that for small Pm in simulations the magnetic Reynolds 
number of the system Rm = Pm Re is in the vicinity 
of the critical magnetic Reynolds number Rm^rit. No¬ 
tice, that Rmcrit ~ 441 for -d = 0.45. At Pm = 0.2 and 
for Re = 1600, the magnetic Reynolds number would 
be Rm = 320 and already be below Rmcrit- A direct 
comparison of our analytical model with the numerical 
simulations is thus difficult in the regime of small Pm. 


IV. CONCLUSIONS 

We present the first detailed scale-dependent model 
for the saturation of the turbulent dynamo. While it was 
known before that this mechanism can efficiently con¬ 
vert turbulent kinetic energy into magnetic energy dur¬ 
ing a phase of exponential growth, the final amplifica¬ 
tion stages and in particular the saturation levels were 
less clear. In this paper we present a phenomenologi¬ 
cal model for the dynamo saturation which generalizes a 
single-scale ansatz suggested by Subramanian [ST] that 
includes the back reaction of the Lorentz force on the 
magnetic diffusivity. We assume that the term describing 
this back reaction is inverse proportional to the turnover 
time of the turbulent eddies and proportional to the mag¬ 
netic energy density. By comparing the dynamo growth 
rate on the scale of fastest field amplification with the 
dissipation rate on the resistive scale we find an expres¬ 
sion for the magnetic energy density. The latter peaks on 
the effective resistive spatial scale which increases in the 
non-linear dynamo phase shifting the magnetic energy to 
smaller wavenumbers. 

In our semi-analytical model we distinguish the cases of 
large and small magnetic Prandtl numbers Pm which are 
fundamentally different due to the different spectral po¬ 
sitions of the resistive and the viscous scale, and ty. 
The two cases are discussed separately in sections [ml 
and |III C] We find a strong dependence of the saturation 
levels on the slope of the initial turbulence spectrum tt. 
The ratio of magnetic energy rrisat over kinetic energy at 
saturation lies between the extreme cases of rrisat/^o and 
Wsat/(eo ~ rrisat), where eo is the initial kinetic energy. 
We expect that the dynamo does not effect the kinetic 
turbulence energy drastically and that eg is maximally 
reduced by the saturation magnetic energy. Our main 
results are listed below: 

• The saturation energy levels are summarized in ta¬ 
ble H For the case of Pm ^ 1 we find values be¬ 
tween 43.8% and 1.3%, while the values for Pm 1 
are with 2.43% to 0.135% considerably lower. The 
saturation level as a function of the turbulence 
spectrum i9 is plotted in figure]^ 

• The evolution of the magnetic spectra is shown in 
figure for Pm ^ 1 and in figure [^ for Pm ^ 
1. By different representative spectra we illustrate 
how the magnetic energy is shifted towards smaller 
wavenumbers during the non-linear dynamo phase 
until saturation is reached. 

• For Kolmogorov turbulence (i9 = 1/3) the peak of 
the magnetic spectrum fc* lies at 101 for Pm ^ 1 
and at 2 for Pm ^1. If the turbulence is a 
highly compressible Burgers type (i? = 1/2) the 
peak scale is at 588 ki, for Pm 1 and at 6 kL for 
Pm ^ 1. 

Our model has some caveats concerning the exact shape 
of the energy spectra. The spectral distribution of mag- 
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netic energy (20) has explicitly been derived for the kine¬ 
matic dynamo phase, when there is no back reaction from 
the magnetic field on the velocity field. With the Lorentz 
force acting on the gas, the spectral shape could in prin¬ 
ciple change. Here we make the simplest possible ansatz 
that the magnetic spectrum remains similar also in the 
non-linear dynamo phase. Another assumption is that 
the slope of the spectrum towards small wavenumbers is 
^3/2 regardless of the slope of the turbulence spectrum. 
As the spectrum has been derived under the assump¬ 
tion of divergence-free turbulence, i.e. Kolmogorov tur¬ 
bulence, it is not clear if this so-called Kazantsev slope re¬ 
mains for compressive turbulence. While there is first nu¬ 
merical evidence that the Kazantsev slope may also hold 
in the compressive turbulence [23], the detailed deriva¬ 
tion of the magnetic energy spectrum for a dynamo in 
this turbulence regime remains an interesting study for 
the future. 

The results from this work can be compared to the satu¬ 
ration model suggested by Schekochihin et al. [53] . These 
authors explain the saturation mechanism of the small- 
scale dynamo in incompressible and non-helical turbu¬ 
lence as a result of anisotropization of the velocity field 
with respect to the magnetic field. While they can re¬ 
produce the spectra from numerical low-Re simulations 
very well, they predict a scaling of the magnetic energy 
at saturation with Re~^^^ leading to tiny values for the 
extremely large Reynolds numbers in astrophysical envi¬ 
ronments. The saturation levels found with our model do 
not scale with the hydrodynamical Reynolds number and 
are moreover independent of the magnetic Prandtl num¬ 
ber in the regimes Pm <C 1 and Pm ^ 1. The expected 
field strength from these saturation levels therefore need 
to be considered in models of the non-linear closures of 
MHD turbulence, and require to consider MHD effects 
even in an initially weakly magnetized gas m- 
An interesting alternative approach for studying the tur¬ 
bulent dynamo are so-called shell models. With an origin 
in hydrodynamic turbulence theory these models have 
been extended to MHD [^. Here the Navier-Stokes 
equation and the induction equation are approximated 
by dividing them into a finite number of shells in Fourier 
space |S^. This allows for solving the MHD equations for 
realistic Reynolds numbers and in particular also in the 
regimes Pm <?; 1 and Pm ^ 1. The exponential growth 
in the kinematic dynamo phase as well as the shift of the 
peak of the magnetic energy spectrum in the non-linear 
phase are well reproduced by shell models [50] ■ The 
magnetic energy spectra in the non-linear phase, how¬ 
ever, seem to deviate from the kinematic shape given in 


equation (20) at high wavenumbers. We also stress that 
Stepanov and Plunian [60] report an excess of magnetic 
energy over kinetic energy not only for high Pm which is 
in agreement to our findings, but within a certain part of 
the inertial range also for low Pm.S 
Section [Thd] of this paper is dedicated to a comparison 
of the predictions from our theoretical model with the 
results found in high-resolution numerical MHD simu¬ 
lations. The values found in simulations by Federrath 
et al. [35| are roughly 40% for Kolmogorov-type turbu¬ 
lence and roughly 1% for Burgers-type turbulence. These 
results are very similar to the predictions of our model 
which are indicated by the blue bands in figure [^ The 
saturation levels as a function of the magnetic Prandtl 
number reported in Federrath et al. |29j are presented 
in figure [^ We find good agreement with our model 
in the large Prandtl number regime for an initial turbu¬ 
lence spectrum with a slope of = 0.45. A quantitative 
comparison with the simulations at Pm ^ 1 is not possi¬ 
ble as the magnetic Reynolds number approaches Rmcrit 
leading to a rapid decrease of the saturation levels with 
decreasing Pm. 

With our phenomenological model we are able to describe 
the saturation of the turbulent dynamo and predict the 
fraction of kinetic energy that is converted into magnetic 
energy. In addition, we can model the magnetic spec¬ 
trum to analyse how the energy is distributed over spa¬ 
tial scales. While numerical simulations are restricted to 
a regime of relatively small Reynolds and Prandtl num¬ 
bers, our analytical saturation model provides for the 
first time estimates for MHD parameters typical for as¬ 
trophysical environments. 
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